Spin dynamics in a one-dimensional ferromagnetic Bose gas 
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We investigate the propagation of spin excitations in a one-dimensional (ID) ferromagnetic Bose 
gas. While the spectrum of longitudinal spin waves in this system is sound-like, the dispersion of 
transverse spin excitations is quadratic making a direct application of the Luttinger Liquid (LL) 
theory impossible. By using a combination of different analytic methods we derive the large time 
asymptotic behavior of the spin-spin dynamical correlation function for strong interparticle repul- 
sion. The result has an unusual structure associated with a crossover from the regime of trapped 
spin wave to an open regime and does not have analogues in known low-energy universality classes 
of quantum ID systems. 
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Quantum interacting onc-dimcnsional (ID) systems 
have for many decades been a subject of a ceaseless inter- 
est of both theorists and experimentahsts. This is mainly 
because of the unique role of quantum fluctuations, which 
are so strong in ID that even for weakly interacting sys- 
tems the intuition based on the free-particle picture and 
the mean field theory fails and the effects of strong cor- 
relations become important Such effects have been 
encountered in many experiments dealing with ID con- 
ductors, like, e.g., organic salts, quantum wires or carbon 
nanotubes, in which the constituent particles, electrons, 
are spin 1/2 fermions. 

Recent advances in the creation and manipulation of 
ultracold atomic gases Q opened an access to a new 
class of ID systems where the constituent particles obey 
bosonic statistics [1, 0, @, @, 0] and have a variable num- 
ber of internal ("spin" ) states In the spinless 
case the theory predicts an equivalence between the Bose 
and Fermi systems: both are described by the Luttinger 
Liquid (LL) theory at low energies In the presence 
of spin the situation is more complex. The necessary 
condition for the applicability of the LL theory is the 
linearity of the dispersion of low-lying elementary exci- 
tations, e{p) ~ \p\. This is usually the case for fermions, 
which have a natural tendency to anti-ferromagnetic or- 
dering [2^. While under special conditions linear dis- 
persion relations, and, consequently, the LL physics can 
be encountered in a multi-component Bose system 
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there exists a broad range of Hamiltonians, in particular 
those with spin-independent interactions, whose ground 
state is ferromagnetic [l^- In the latter case the soft- 
est low-lying excitation is the magnon with a quadratic 
dispersion relation 



e{p) ~ p^/2to. 



(1) 



where m» is an effective mass. This makes the straight- 
forward application of the LL theory impossible and 
poses a fundamental question of finding an alternative 
theory describing the dynamics of the low-energy excita- 



tions in a ID boson ferromagnet. 

In this Letter we tackle this issue and compute the 
long-distance properties of two-point correlation func- 
tions of local spins at zero temperature. Focusing on 
the region of strong interparticle repulsion, we show that 
the dynamical properties of spin excitations in a ID fer- 
romagnetic Bose gas are neither those of a localized fer- 
romagnet nor of a Luttinger liquid, pointing at the exis- 
tence of a new low-energy universality class. Our main 
results are presented in Eqs. (fT3|) and p^ . We also dis- 
cuss the connection of our work with the problems of 
a moving impurity in a LL, dynamics of a hole in the 
Hubbard-Mott insulator, and quantum mechanics in a 
dissipative environment. Finally, we describe recent ex- 
perimental realizations [1, 0, [13] of quasi-lD Bose gases 
with spin. 

For simplicity of presentation the derivations are car- 
ried out for two-component bosons; the generalization to 
higher spins is straightforward [l3|- We assume that the 
interaction between the particles is spin-independent and 
the model Hamiltonian has the form 



H 



^ 2m 



gS{a 



,) + U{x,^x,)]+hS,. (2) 



Here m is the mass of a boson, h is the external mag- 
netic field, and Sz is the z component of the total spin. 
We are interested in the limit of infinite number of par- 
ticles, N — > 00, and of infinite system size, L — > 00, at 
a fixed particle density, po = N/L. Although for cold 
atoms the interaction potential is well approximated by 
a delta- function, we allow for a more general interaction 
gS{x) + U{x), where U{x) is some smooth function. The 
strength of the short-range repulsion is characterized by a 
dimensionless coupling constant 7 = mg/h?po. For U = 
the Hamiltonian ^ can be diagonalized by Bethe Ansatz 
(BA) [1^ , providing us with a valuable source of intuition 
about the low-energy dynamics studied here. 

The global spin operator S = {Sx, Sy, Sz) can be rep- 
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FIG. 1: A ID array of particles carrying spin. The propa- 
gation of longitudinal (a) and transverse (b) spin waves over 
the fully polarized state | -^1) is depicted. The state obtained 
by the action of the spin lowering operator s_(0, 0) onto | -[]■) 
is schematically illustrated in (c). 



resented as 



S = 



/ dxs{x), s(a;) = -^, 
^0 ^ ^-^1 



d{x-Xj), (3) 



where cr^-'^ = {ux, cFy, ctz) is the vector composed of the 
three Pauli matrices acting non-trivially on the spin in- 
dices of the j-th particle. Spin-ladder operators s±{x) = 
Sx{x) ± isy{x) flip the z component of a local spin. For 
/i > the Hamiltonian ^ has a non-degenerate ground 
state, I fl"), which is fully polarized along the z-axis, 
s+(x)| ff) = 0. The degeneracy appears at /i = and 
is discussed in detail in Ref . [IJI . We investigate the dy- 
namics of excitations over the state | f|~). For simplicity 
of the presentation, our results are given for ft, = 0. In 
the case of ft 7^ the r.h.s. of Eqs. ([6]) and should 
be multiplied by the oscillating factor e'*'*. 

Our main object of interest is the response of the sys- 
tem to weak perturbations of local magnetization, en- 
coded into two correlation functions: 

G||(.T,t) = (^|s,(a:,t)s,(0,0)|^) (4) 
Gi(x,t) = (^|s+(x,t)s_(0,0)|^) (5) 

describing the longitudinal. Fig. [Ija), and transverse. 
Fig. [TJb), spin dynamics, respectively. The longitudi- 
nal dynamics over the state | f|-) coincides with that of 
density fluctuations in the ID spinless Bose gas. Their 
dispersion is linear in the low-energy limit, therefore the 
X, t — > cx) asymptotics of Eq. (jj]) can be calculated within 
the LL theory |H, and decays as a power-law [27| . 

In contrast to the longitudinal dynamics, transverse 
spin waves are not sound waves, their dispersion is 
quadratic at low momenta, Eq. ([T]). The effective mass 
m* in Eq. ([1]) increases with increasing 7. In the BA solv- 
able case m* can be calculated exactly 15| and shows 
a linear divergence with 7 in the limit of strong repul- 
sion: ra^jm ~ 87/271^ as 7 — > 00. This divergency was 
pointed out in Ref. [l^l as a signature of a slow dynamics 
of transverse spin waves at large 7. We show here that 
in the strong coupling regime the effects are even more 
dramatic. To see this, we start with the description of 
the long-wavelength limit of Eq. ([5]) at 7 = 00. 



Qualitatively, the propagation of the transverse spin 
excitation at 7 = cx) can be described as follows: the 
operator s_(0,0) flips a spin of a given particle, shown 
in Fig. [ijc) with a down arrow. Because of the infinite 
repulsion, the spin-down particle cannot exchange its po- 
sition with its neighbors, so it is trapped inside an inter- 
val of average length 1p^^ ■ However, since the spin- up 
particles are mobile, the size and the position of this in- 
terval can fluctuate, making it possible for the spin-down 
particle to be detected at a distance larger than p^^ from 
its initial position. Such fluctuations require a simultane- 
ous displacement of a large number of spin-up particles, 
thus the probability to observe the spin-down particle at 
a large distance from the point where it was created must 
be small. The correlation function ([5]) should thus decay 
very rapidly with distance. 

To quantify the above argument we obtain the ana- 
lytic expression for the long-wavelength asymptotics of 
Eq. (|5l). The calculations are performed using, as in 
Ref. [l6|, a combination of bosonization with a first- 
quantized path integral. For the BA solvable case they 
can be underpinned by the asymptotic analysis of the 
determinant representation of the correlation functions 



along the lines of Refs. [17|, |18[. We give only the result 



here, the details will be presented elsewhere [l3|. The 
asymptotic form of Gx(a;,i) is given for t ^ by: 



Gi(.T,t) 



1 



: CXp 



^_{J^)Qxf_ 

' K2\n{t/tv) 



(6) 



It is similar to the diffusion propagator except that the 
mean square deviation grows logarithmically with time. 
This logarithmic diffusion is the mathematical manifes- 
tation of the spin trapping effect. The parameter con- 
trols the smallest time scale in the problem: 



h 



2m 



(7) 



The physical meaning of tp can be understood from 
the fact that the spectrum of the (5-interacting spinless 
bosons at 7 = cx) is the same as that of free spinless 
fermions. Therefore, Ep in Eq. ([7]) plays the role of a 
Fermi energy, and tp is the time scale for the longitu- 
dinal spin fluctuations. The dimensionless parameter K 
in Eq. ^ is the Luttinger parameter, which can be cal- 
culated from the thermodynamic properties of the sys- 
tem Note that for ?7 = one has K > I and K 1 
only at 7 — !■ 00. 

Eq. ([6]) is obtained for 7 = cxo. For large but finite 7 
there is a small probability that spin-up and spin-down 
particles exchange their positions, allowing the spin ex- 
citation to escape from the trap. We estimate the es- 
cape time by replacing the fluctuating gas of spin-up 
particles by a static Kronig-Penney lattice with period 
Pq ^ and get jtp. Thus, for large 7 the escape time 
is parametrically large and there exists a broad window 
tp ^ t t^, where G±{x, t) has the asymptotics form 
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For t >^ tif, there is a crossover to another, "open", 
regime, which we investigate next. A rather general 
result can be obtained from a minimal set of assump- 
tions on the analytic properties of the spectral function 
A{k,uj). By definition: 



G^{x,t) 



dk 

£ 

27r 



— ( 

2-K 



*A{k,uj) 



with 



A{k,oj) = Y,5{nw - E,{k))\{v, k\s-{k)\ ^)|^ 



(8) 



(9) 



Here H\v,k) — E^{k)\v^k) and v enumerates all the 
states with a given momentum hk. Among these states 
there is a state of a minimal energy e(fc), that is e{k) = 
min^ E^{k). Our first assumption, supported by the cal- 
culations for the BA solvable case [13|, where e{k) can 
be found explicitly, and variational considerations in the 
spirit of the Feynmann single-mode approximation is that 
at small p = hk the function £{k) has the form ([T]). 
Therefore, £{k) defines the threshold frequency for the 
spectral function: A{k,uj) = for ftw < e{k). Above the 
threshold a continuum of states contributes to A{k^Ljj). 
They contain one magnon (spin-flip) and multiple plas- 
mon (density) excitations close to zero momentum. The 
large t asymptotics of Eq. ^ is dominated by the scale- 
free part of the spectral function at the threshold, whose 
most generic form is 



A{k,u) ~ c{k)[nw - e{k)]^^''\ hio > e{k), 



(10) 



where A(fc) and c{k) are some functions of k. 

Our second assumption is that both c(fc) and A(fc) are 
analytic functions of momentum in the vicinity of fc = 0. 
Taking into account the inversion symmetry, G_\_{x^t) = 
G±{—x,t), which imphes A{k,uj) = A{—k,uj), the gen- 
eral form of c(fc) and A(fc) should be 



A(fc) = a-l + /3fc2 
c(fc) = Co + cifc^ + • • 



(11) 
(12) 



where a, (3, cq and ci are some model-dependent co- 
efficients. We substitute Eqs. ([TO]), ^ and ^ into 
Eq. ([H]), and for the function e(fc) entering Eq. pU]) we 
use formula ([T]) with p = hk. The saddle point analysis 
of the resulting expression leads to our main result: 



Gj.(x,t) ~t-" 



, t \ ith 
f3ln{ — ]+ 



-1/2 



X exp 



2th-Ail3m^ ln(t/tF) 



(13) 



Eq. p3p contains two parameters, a and /3, which depend 
on both 7 and U. This dependence cannot be extracted 
from the scaling and analyticity assumptions, Eqs. ([TT 



and (|lip . However, we can find a and /3 in the trapped 
regime, 7 = cxd, by letting m* = 00 in Eq. ([13]) and 
comparing the resulting expression with Eq. ([6]) . We get 



0, 



K 



2(7rpo)^ 



(14) 



at 7 = 00. In addition, for [/ = we could treat the case 
of large but finite 7 by BA and showed that Eq. pi]) 
remains valid 13|. Whether this is an artefact of inte- 



grability or a generic property of the model is an im- 
portant open question. The latter would mean that the 
asymptotic behavior of the correlation functions in the 
ferromagnetic case is completely determined by thermo- 
dynamic properties of the system, like in the LL theory. 

If m, < 00 equation coincides with ^ for suffi- 
ciently small times, 2th ^ 4/?TOh. ln(t/tF). At a time 
when this relation becomes an equality a crossover to an 
open regime occurs For f/ = the dependence of 
on 7 can be found explicitly: 



(15) 



We stress that this estimate is valid if 7 is large enough to 
ensure 3> tp- This confirms and completes our naive 
estimate of given in the paragraph below Eq. ([6]). 

In the open regime, t >tt, the function factorizes 
into a product of the transverse spin correlation func- 
tion of the localized Heisenberg ferromagnet, G^, ex- 
hibiting exchange-induced oscillations, and a rapidly de- 
caying factor, which we attribute to the excitation of the 
charge degrees of freedom: 



H 

± J 



m 



2K- 



(16) 



The strong suppression of G± for x > l{t) ensures the 
absence of the light-cone singularities at x = vt, where 
1; is the sound velocity for the charge excitations, v = 
hnpQ/m for 7 = 00. This is illustrated in Fig. [21 

We now discuss the relation of our work to several 
previously studied problems. The threshold singular- 
ity, Eq. ([TO)) , of A{k,u)) is identical to the Fermi edge 
singularity in a ID quantum fluid with a mobile impu- 
rity [3, [i^, [2lj in a special case of equal impurity and 
host particles masses. Our results, therefore, have inter- 
esting implications for the latter problem. In particular, 
the relation of A(fc) to the parameters of the model has 
only been understood at fc = 1^, 2^ . The finite fc case 
was addressed in Ref. [2l[ by using an effective field the- 
ory. We, however, found that A(fc) calculated from the 
Ref. [21} disagrees with our Eqs. ([TT|) and ([T4)) . In partic- 
ular, it does not reproduce Eq. ([6|), which we obtain by 
two independent methods [2^. There is also an analogy 
with the problem of a mobile hole in a ID Hubbard-Mott 
insulator (2^ . There the parabolic branch of spectrum 
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FIG. 2; (color online) Shown is the intensity plot of 
t'^Re G± (x, t) in the x,t plane at 7 = 100. The spin wave 
is strongly suppressed at a distance x ~ (-{t), Eq. (|16p . long 
before it reaches the light cone, x — vt, shown as dashed 
green lines. The onset of the spin precession is seen in the 
space oscillations of t"IieG± developing above the crossover 
time t* ~ IStp, shown by the white dotted line. 



is due to holon excitations, while the gapless mode with 
a linear dispersion corresponds to spinons. The spectral 
function given by Eq. (3) of Rcf. ^22] shows the same crit- 
ical properties as A{k,uj), Eqs. (fTU)) and pT|) . We thus 
hope that the results presented in this Letter can be ex- 
tended to a broader class of models. Finally, by bosoniz- 
ing the gas of spin-up particles one maps the Hamilto- 
nian ^ onto that of a quantum particle non-linearly cou- 
pled to a harmonic environment. A linearized version of 
this problem was studied in Refs. [23, 24|, and the com- 
parison of the long-time asymptotics remains an open 
question. 

Experiments on quasi-lD Bose gases with internal de- 
grees of freedom started recently [^, [l^ . In the experi- 
ment Q a pseudospin 1/2 system was created by loading 
®^Rb atoms into a highly elongated magnetic trap and 
selecting two hyperfine states, IF = l,mF = —1) and 
|F = 2,mF = 1), while in [3, [l^ the true spin 1 sys- 
tem was created by loading *^Rb atoms into a highly 
elongated optical trap and selecting the states compos- 
ing F = I spin triplet. In both Q and [1, [l^l long-lived 
fully polarized states were achieved and a possibility to 
excite and image (pseudo)spin waves in real space and 
time was demonstrated. Although in the above exper- 
iments the systems were not truly ID (about 50 bands 
of transverse quantization were populated) a possibility 
to further reduce the number of occupied bands is sug- 
gested by the successes in the creation of ID spinless 
systems 0, 0, H, @, 01 : ttius paving the way to the inves- 
tigations of the ferromagnetism in ID Bose gases. 
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